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ABSTRACT 


The  purpose  of  this  study  was  to  indentify  preferred  areas  of  synoptic 
scale  disturbance  activity  through  the'/use  of  the  objective  analysis 
technique  called  complex  principal  component  analysis.  The  analysis 
technique  was  applied  to  a  sea  level  pressure  data  set  consisting  of 
twice  daily  observations  (00  and  12  GMT)  during  the  months  of  January, 
February,  and  March,  for  the  years  1973-1976  and  1979- 1 982 .  A  Ttre  months 
of  January  through  March -were  chosen  to  encompass  the  time  frame  of  the 
Genesis  of  Atlantic  Lows  Experiment  1986.  The  geographical  area  studied 
extended  approximately  from  30-120  degrees  west  and  from  25-60  degrees 
north. 

The  complex  principal  component  analysis  technique  allows  the 
detection  of  propagating  features  in  the  pressure  data  set.  The  original 
data  set  is  augmented  with  its  own  Hilbert  transform,  which  is  nearly 
equivalent  in  the  time  domain  to  the  quadrature  spectrum  in  the 
frequency  domain.  ,  . 

The  variance  in  the  pressure  field  is  felt  to  be  a  direct 
reflection  of  disturbance  activity.  Furthermore,  when  the  data  are 
filtered  to  a  "synoptic  scale"  frequency  window,  the  variance  in  the 
pressure  field  corresponds  closely  with  cyclone  and  anticyclone  storm 
tracks  over  a  time  period  of  two  to  ten  days. 

The  results  of-the  analysis  produced  four  components  that  could  be 
given  a  preliminary  interpretation  from  a  physical  sense.  The  components 
(patterns)  resembled  areas  of  activity  as  discovered  in  previous  studies 
using  other  subjective  analyjiS  techniques. 
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1.  INTRODUCTION 


Considerable  time  and  effort  have  been  devoted  to  the  study  of 
synoptic  cyclones  and  anticyclones  and  their  preferred  areas  of  genesis 
and  movement.  Synoptic  scale  cyclones  (lows)  and  anticyclones  (highs) 
have  wavelengths  of  1000  to  2500  km,  as  seen  in  the  lower  troposphere. 
They  are  recognized  on  synoptic  surface  weather  charts  as  migratory 
low-  and  high-pressure  systems,  respectively.  This  keen  interest  in 
synoptic  scale  disturbances  is  understandable  since  the  daily  weather 
we  experience  is  closely  related  to  the  cyclones  and  anticyclones  that 
form,  move  across,  and  dissipate  over  the  earth's  surface.  These 
disturbances  are  features  of  the  general  circulation  and,  as  such,  have 
preferred  areas  of  activity  as  controlled  by  the  general  circulation. 

As  indicated  by  Barber  et  aK  (1984)  there  have  been  many  studies 
of  the  spatial  distribution  of  cyclone  activity  (anticyclone  studies 
are  fewer  in  number),  and  a  few  studies  that  examined  the  temporal 
variations  in  synoptic  scale  disturbance  activity.  Most  of  these 
studies  are  based  on  counts  of  individual  disturbances  on  daily  weather 
maps  or  on  published  disturbance  track  maps.  These  subjective  counting 
techniques  have  some  inherent  problems  associated  with  them,  including 
changing  methods  of  map  analysis  and  interpretation,  changing  data 
density  over  time,  and  unequa"  grid  areas  with  changing  latitude. 

To  overcome  these  problems  an  objectively  determined 
meteorological  variable  is  needed  which  is  accurately  measured  and 
whose  variance  is  related  to  the  synoptic  scale  disturbance  activity. 
Barber  et  ajL  (1984)  chose  the  station  pressure  as  the  variable 


2 


required.  It  is  accurately  measured,  and  studies  by  Klein  (1951a, 
1951b),  Blackman  et  aK  (1977),  Hartmann  (1974),  and  Barber  et  al . 
(1984)  strongly  suggest  that  there  is  a  spatial  correspondence 
between  synoptic  scale  disturbances  and  the  variance  in  station 
or  sea  level  pressure,  especially  if  the  pressure  data  are  filtered  to 
look  at  an  appropriate  "synoptic  scale"  frequency  window  with  a  period 
of  approximately  2-10  days. 

Barber  et  aK  (1984)  applied  spectral  analysis  techniques  to  a 
limited  data  set  (3  data  stations)  to  test  the  usefulness  of  the 
station  pressure  as  an  indicator  of  disturbance  activity  and  found 
further  support  for  the  relationship  between  synoptic  scale 
disturbances  and  the  variance  in  the  station  pressure  field.  The 
recently  documented  analysis  technique  called  complex  principal 
component  analysis  (Barnett,  1983;  Horel ,  1984)  appears  to  be  a 
promising  technique  to  apply  to  a  geographically  larger  data  set  of 
pressure  values  to  look  for  relationships  between  the  variance  in  the 
pressure  field  and  synoptic  scale  disturbance  activity. 

Complex  principal  component  analysis  (CPCA)  in  the  time  domain  is 
designed  to  take  a  set  of  linearly  interdependent  variables,  which  are 
functions  of  space  and  time,  and  transform  (rotate)  the  original 
configuration  of  observations  to  a  new  orientation  having  the  same 
number  of  dimensions  (variables)  but  where  the  new  orientation  consists 
of  a  set  of  mutually  orthogonal  (independent)  dimensions.  These  new 
dimensions  (components)  are  linear  composites  of  the  original  variables 
and  account  for  the  variation  in  the  original  variables.  The 
components  accounting  for  the  largest  amounts  of  variance  in  the 


original  data  can  be  used  to  study/interpret  the  physical  processes 
taking  place  in  the  original  data  set.  Since  the  components  are 
orthogonal,  they  can  be  combined  in  various  ways  or  studied  separately. 
Sometimes  a  large  percentage  of  the  variation  in  the  original  variables 
can  be  accounted  for  by  the  first  few  largest  components.  This  is  the 
goal  of  CPCA--to  reduce  the  number  of  dimensions  needed  to  "explain"  a 
data  set  using  new  dimensions  that  are  linearly  independent. 

CPCA  allows  the  efficient  detection  of  propagating  features  in  a 
geophysical  data  field  that  are  derived  from  the  interaction  between 
traveling  and  standing  waves  of  different  spatial  scales  and  temporal 
frequencies  (Horel,  1984).  CPCA  should  allow  the  separation  of  the 
preferred  patterns  (areas)  of  synoptic  scale  disturbances  from  the 
overall  variance  in  the  sea  level  pressure  field  since  the  individual 
disturbances  are  guided  by  the  standing  and  traveling  waves  in  the 
general  circulation  patterns. 

The  sea  level  pressure  was  chosen  as  the  variable  to  use  in  this 
study  versus  the  station  pressure  used  by  Barber  et  aK  (1984).  By 
using  the  sea  level  pressure  the  elevation  of  the  observation  stations 
is  eliminated  as  a  variable  in  the  computations.  Standard 
observational  practice  is  to  reduce  the  observed  station  surface 
pressure  to  a  value  that  would  exist  at  a  point  at  sea  level  directly 
below  the  station,  if  the  temperature  at  sea  level  is  assumed  to  be  the 
same  as  the  temperature  at  the  higher  elevation  station.  The  actual 
method  used  to  reduce  the  pressure  uses  the  preceding  12-hour  mean 
temperature,  not  the  instantaneous  temperature  at  the  time  of  the 
station  pressure  observation.  This  introduces  some  error  into  the 
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pressure  observations,  but  not  enough  to  outweigh  the  advantages  of 
eliminating  the  elevation  as  a  variable. 

The  objectives  of  this  study  are  to  1)  apply  the  objective 
analysis  technique,  CPCA,  to  a  sea  level  pressure  data  set  and  try  to 
detect  the  preferred  areas  and  movement  of  synoptic  scale  disturbances; 
2)  document  the  areas  and  movement  discovered  and  show  their  agreement 
or  disagreement  with  previous  study  results;  3)  and  explain  the 
limitations  and  problems  encountered  during  this  study  effort  and  how 
they  affected  the  results  of  this  analysis. 
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2.  PRELIMINARY  DATA  PROCESSING 

The  data  used  in  this  study  are  the  final  analyses  of  the  sea 
level  pressure  field,  for  the  Northern  Hemisphere,  as  analyzed  by  the 
U.S.  Navy  using  the  objective  analysis  scheme  called  "Fields  by 
Information  Blending"  ( Hoi  1  and  Mendenhall,  1971).  The  data  set  was 
obtained  from  the  data  archives  of  the  National  Center  for  Atmospheric 
Research  (NCAR).  There  are  two  maps  per  day,  at  00  and  12  GMT,  for  the 
time  periods  from  1  January  to  31  March,  for  the  years  1973  to  1976  and 
1979  to  1982. 

The  NCAR  data  set  contains  sea  level  pressure  values  (hereafter 
referred  to  as  SLP)  interpolated  to  an  octagonal  grid  system  which  is  a 
square  grid  overlaid  on  an  octagonal  area  of  a  polar  stereographic 
projection  of  the  earth’s  surface.  The  grid  points  are  approximately 
324  km  apart,  depending  on  the  latitude  (for  further  details,  see 
Jenne,  1970,  1975).  See  figure  2 a  for  the  geographical  area  covered  by 
the  data  set  used  in  this  study. 

The  time  period  of  1  January  to  31  March  was  selected  to  encompass 
the  time  frame  of  the  Genesis  of  Atlantic  Lows  Experiment  (GALE)  1986. 
The  preferred  areas  of  synoptic  disturbance  activity  and  their  movement 
during  the  winter  months  were  of  interest  for  GALE  purposes.  Of 
particular  interest  are  the  cyclones  that  develop  off  the  Atlantic 
coast  and  intensify  rapidly  as  they  move  northward.  The  data  set  years 
of  1973-1976  and  1979-1982  were  selected  based  on  the  fact  they  had  the 
least  number  of  observational  maps  missing  (approximately  1  percent), 
and  no  more  than  one  map  missing  consecutively.  This  allowed  the 


replacement  of  missing  maps  with  the  values  obtained  from  a  linear 
interpolation  from  neighboring  values  in  time. 

The  mean  SLP  at  each  individual  grid  point,  for  January  through 
March  for  each  year,  was  calculated  and  subtracted  from  the  individual 
SLP  observations  in  the  respective  yearly  time  series.  The  time  series 
were  adjusted  slightly  to  compensate  for  leap  years  in  1976  and  1980. 
The  SLP  deviations  were  then  Fourier  transformed  and  the  Fourier 
coefficients  for  periods  greater  than  10  days  or  less  than  2  days  were 
set  equal  to  zero.  The  inverse  transform  resulted  in  a  band-passed 
filtered  time  series  with  information  only  in  the  desired  time  interval 
of  synoptic  scale  events. 

The  period  described  above  was  based  on  previous  work  which  showed 
disturbances  in  this  time  interval  "exhibited  variance  and  covariance 
patterns  characterized  by  elongated  maxima  suggestive  of  storm  tracks 
and  an  association  with  baroclinic  waves"  (Blackmon  et^  aK ,  1984a),  and 
that  these  disturbances  contribute  the  most  to  the  spectrum  of  the 
variance  when  studying  the  time  period  from  .25  to  18.75  days  (Barber 
et  aK ,  1984).  Other  studies  by  Blackmon  (1976),  Blackmon  and  White 
(1982),  Blackmon  et  al_.  (1977,  1979,  1984a,  1984b),  White  (1980),  and 
Hartmann  (1974)  all  suggest  a  spatial  and  temporal  correspondence 
between  cyclones/  cyclone  tracks,  and  variance  of  SLP  and  other 
parameters  at  upper  levels  in  the  atmosphere  (e.g.,  850  mb  temperature 
deviations;  700,  500,  300  mb  height  deviations).  Care  must  be 
exercised  as  regions  of  high  variance  may  also  correspond  to  tracks  of 
anticyclones  as  well  as  cyclones  (Klein,  1951b),  and  this  was  found  to 
be  the  case  as  discussed  later. 


3.  BACKGROUND  THEORY  AND  METHODS 


The  development,  in  the  meteorological  literature,  of  CPCA  in  the 
time  domain  can  be  traced  back  to  Rasmusson  et  al_.  (1981),  while  CPCA 
in  the  frequency  domain  appeared  first  in  Wallace  and  Dickinson  (1972) 
In  this  section  a  brief  summary  of  relevant  theory  will  be  presented 
along  with  the  methods  used  in  this  study. 


a.  decomposition  theory 

The  work  of  Rasmusson  et  aK  (1981)  was  based  on  the  Singular 
Value  Decomposition  (SVD)  theorem.  The  theorem  supports  the 
relationship  between  an  arbitrary  real  nxm  (n  <  m)  matrix  A  =  {A.  .} 

*  J 

and  a  summation: 


A  -  ;  +  x  x<2M2)  + 

Aij  lXi  \j  A2Xi  Yj 


+  A  X^M"^ 
n  l  j 


with  >i  >  *2  >  •••  *n  >  0,  where  the  X  * s  are  referred  to  as  singular 
values  and  the  {x(£^}  ^  n  and  m  are  known  as  the  nth  left 
and  right  singular  vectors,  respectively. 

The  theorem  can  be  expressed  in  matrix  form  as 


A  =  X  a  Y'  (2) 

nxm  nxr  rxr  rxm  ’  v  ' 

where  X  and  Y  are  each  orthonormal  sections  and  a  is  a  diagonal  matrix 
of  ordered  positive  singular  values.  The  prime  denotes  the  transpose. 
The  rank  of  A  is  given  by  r,  the  number  of  positive  values.  The 
columns  of  X  and  Y  are  the  left  and  right  singular  vectors, 
re^pecti vely,  of  A.  The  left  singular  vectors  (columns  of  X)  are  the 
eigenvectors  of  the  product  AA' ,  and  the  right  singular  vectors 
(columns  of  Y)  are  the  eigenvectors  of  the  product  A'A.  Thus,  the 


eigenvectors  are  derived  from  A.  Each  of  these  product  matrices  has 
eigenvalues  given  by  the  respective  singular  value  squared.  That  is, 
the  matrix  Arxr  is  a  diagonal  matrix  whose  main  diagonal  entries  are 
the  square  roots  of  the  positive  common  eigenvalues  of  AA '  and  A'A 
(Green,  1978). 

If  the  original  matrix.  A,  which  was  n  x  m,  had  observations  of  a 
variable  in  time  indexed  by  i  =  1,2,. ..n,  and  grid  point  locations 
indexed  by  j  =  l,2,...m,  then  the  operation  AA'  would  produce  a  n  x  n 
matrix  which  contains  only  information  concerning  the  temporal 
characteristics  of  the  original  matrix  A.  The  operation  A’A  produces  a 
m  x  m  matrix  which  contains  only  information  concerning  the  spatial 
characteristics  of  the  matrix  A.  The  SVD  technique  thus  allows  the 
separation  of  information  describing  the  data  matrix  under  study. 
Additionally,  since  the  A  ^  matrix  contains  singular  values  common  to 
both  the  AA1  and  A'A  matrix  operation  results,  the  most  important  terms 
from  AA'  and  A'A  can  be  selected  to  describe  A  based  on  the  magnitude 
of  the  associated  squared  singular  values. 

The  leading  terms  of  the  expansion,  (1),  constitute  an  efficient 
basis  for  describing  a  given  data  array.  The  terms  can  be  ordered  so 
that  the  terms  with  the  largest  A ' s  are  first  in  the  expansion  since 
the  terms  are  orthogonal  and  the  vectors  have  unit  length.  It  can  be 
shown  that  this  expansion  is  best  in  the  least  squares  sense.  The 
magnitude  of  x  directly  establishes  the  importance  of  each  term  in  the 
expansi on--the  leading  terms  accounting  for  the  largest  portion  of  the 
variance  in  the  given  array,  and  succeeding  terms  accounting  for 
smaller  and  smaller  portions  of  the  variance. 
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The  potential  usefulness  of  this  decomposition  method  is  realized 
when  it  is  applied  to  a  data  set  containing  features  that  are  highly 
correlated  on  broad  scales.  A  few  leading  terms  of  the  decomposition 
(expansion)  will  usually  contain  the  most  meaningful  information 
contained  in  the  matrix,  while  the  remaining  terms  contain  signals 
considered  noise.  It  should  be  easily  seen  that  in  calculating  tlv 
eigenfunction  of  one  product  matrix  Y  (=A‘A),  the  eigenfunction  of  the 
other  product  matrix  X  ( =AA ' )  follows  from  a  simple  matrix 
multiplication  operation  when  using  the  SVD  technique. 

The  SVD  theorem  outlined  above  can  be  altered  to  include  complex 

matrices  also.  All  matrix  transpose  operations,  indicated  by  the  prime 

notation,  would  become  complex  conjugate  transpose  operations. 

Additionally,  the  references  to  in  (1)  would  become  the  complex 

(  e) 

conjugate  of  Yj  '. 

While  the  decomposition  of  a  real  nonsymmetric  matrix  is  usually 
referred  to  as  SVD,  the  decomposition  of  a  real  symmetric  cross 
products  matrix  (i.e.,  correlation  matrix,  covariance  matrix,  etc.)  is 
often  referred  to  as  principal  components  analysis  (PCA). 


b.  principal  components  theorv 


The  SVD  of  a  data  matrix.  A,  can  be  found  with  A  defined  as 


A  *  (Afj  -  Sj)  ,  (3) 

where  i  =  1,2,. ..n  and  j  =  1,2,. ..m,  and  A.  is  the  mean  of  each  column 

J 

of  A.  Define  also  the  covariance  matrix,  C,  as 


A ' A  /  (m-1)  . 


(4) 
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The  SVD  of  a  symmetric  cross  product  matrix  (C  in  this  case)  takes 
the  form 

C  =  UDU'  ,  (5) 

where  U  is  an  orthogonal  matrix  of  eigenvectors  and  D  is  a  diagonal 

matrix  of  C's  eigenvalues.  Thus  it  can  be  seen  that  the  SVD  of  a 

symmetric  product  moment  matrix  provides  (Green,  1978):  1)  a  matrix  of 

2 

eigenvalues  D  in  the  decomposition  of  C  which  is  equal  to  a  /(m-1),  as 
would  be  obtained  from  the  SVD  of  the  matrix  A;  2)  a  matrix  of 
eigenvectors  U  in  the  decomposition  of  C  which  is  exactly  the  same  as  Y 
in  the  SVD  of  A. 

If  the  matrix  product  A'A  has  been  constructed  to  contain  spatial 
information  only,  then  after  computation  of  U  (=Y)  we  can  find  Z  such 
that 

Z  =  AU  =  Xa  ,  (6) 


where  Z  contains  the  temporal  information  originally  contained  in  A. 

We  will  refer  to  the  columns  of  U  (eigenvectors)  as  the  components 
(which  show  spatial  variation)  and  to  the  columns  of  Z  as  the  component 
scores  (which  show  temporal  variation). 

The  components  can  be  viewed  as  sets  of  direction  cosines  that 


the  projection  of  the  original  data  observations  onto  the  component 
axes.  The  scores  indicate  how  important  an  observation  is  in 
interpreting  a  component.  When  the  scores  are  plotted  as  a  time 
series,  the  plot  will  reveal  at  what  time  the  observations  had  the 
greatest  variability.  The  units  of  the  scores  are  the  same  as  the 
units  on  the  original  variables.  The  eigenvalues  contained  in  the 
matrix,  D,  indicate  how  much  of  the  total  variance  of  the  original  data 
set  is  accounted  for  by  each  pair  of  components  and  component  scores. 
The  units  of  the  eigenvalues  are  the  same  as  the  units  on  the  original 
variables,  but  squared. 

Many  standard  meteorological  variables  are  interdependent  in  space 
or  time,  or  both--some  related  in  ways  so  complex  that  physical  cause 
and  effect  links  are  not  clear.  PCA  transforms  a  set  of 
intercorrelated  variables  into  a  new  coordinate  system  in  which  the 
axes  are  linear  combinations  of  the  original  variables  and  are  mutually 
orthogonal  (Hayden  and  Smith,  1982).  The  objectives  of  the  PCA  are  to: 
1)  untangle  complex  patterns  of  intervariable  association  in 
multivariate  data;  2)  isolate  characteristic,  recurrent,  and 
independent  patterns  of  covariance  among  the  variables  into  a  new  set 
of  independent  variables;  and  3)  reduce  the  dimensionality  of  the  set 
of  multivariate  data  and  still  retain  its  most  important  features. 

PCA  can  be  considered  a  closed  system  model  in  that  all  the 
variance  in  the  original  variables  is  being  investigated  (Johnston, 
1978).  The  analysis  technique  does  not  depend  on  any  predetermined 
form,  but  rather  depends  upon  the  interrelationships  within  the  data 
being  analyzed.  The  form  of  the  component's  patterns  is  determined 
empirically  from  the  fundamental  data.  The  variance  of  a  scalar  (or 


vector)  field  can  be  represented  by  comparatively  few  independent 
coefficients  of  the  eigenvectors  when  plotted  spatially,  and  the 
associated  patterns  may  in  some  cases  be  given  a  physical 
interpretation. 

If  the  eigenvectors  are  derived  from  a  covariance  matrix,  their 
spatial  form  depends  directly  on  the  coherence  of  the  departures  from 
normal  over  the  area  of  analysis.  The  eigenvector  patterns  will 
reflect  the  real  spatial  distribution  of  variance  over  the  study  area. 
If  the  covariance  matrix  is  based  on  SLP  data,  the  resulting 
eigenvector  patterns  should  represent  synoptic  patterns.  The  pattern 
represented  by  the  eigenvector  (component)  associated  with  the  largest 
eigenvalue  can  be  thought  of  as  representing  the  mean  pattern  of 
variability  in  the  data  field  (Johnston,  1978). 

PCA  has  a  major  limitation  in  that  it  only  resolves  instantaneous 
temporal  covariance;  thus  an'./ialies  that  move  in  space  and  time  cannot 
be  fully  identified  (Michaelsen,  1982).  It  can  only  portray  variations 
in  time  as  either  in  phase  or  180°  out  of  phase.  So  it  cannot  be  used 
to  represent  features  with  variable  phase  relationships,  such  as 
traveling  waves. 


c.  frequency  domai..  CPCA 

Wallace  and  Dickinson  (1972)  developed  what  has  been  called  CPCA 
in  the  frequency  domain  (hereafter  referred  to  as  FDPCA).  They  were 
trying  to  develop  a  technique  to  avoid  problems  inherent  in 
cross-spectral  analysis  of  data  sets  where  there  are  more  than  one  type 
of  wave  structure  present  in  the  same  frequency  band,  i.e.,  there  is  no 


way  to  determine  how  many  significant  wave  structures  are  present  or 
their  relative  contribution  to  the  variance  spectra. 

Wallace  and  Dickinson's  new  technique  involved  computing  the 
cross-spectrum  matrix  of  the  data  set,  where  the  diagonal  elements  o, 
the  matrix  are  the  power  spectral  estimates  obtained  by  crossing  each 
time  series  with  itself,  and  the  off-diagonal  elements  are  complex,  the 
real  part  being  the  co-spectrum  and  the  imaginary  part  the  quadrature 
spectrum.  The  original  data  set  was  first  band-passed  filtered  to 
eliminate  all  components  outside  the  frequency  band  of  interest.  The 
vector  time  series  are  expressed  as  a  linear  combination  of  the 
eigenvectors  of  the  cross-spectrum  matrix.  The  eigenvectors  are 
applied  to  an  augmented  time  series  involving  the  original  time  series 
and  its  time  derivative. 

Wallace  and  Dickinson  (1972)  showed  that  a  time  series  u  - ( t )  has  a 

J 

spectral  representation  (analogous  to  Fourier  representation), 

oo 

u.(t)  =  Re  /  exp(iwt)  dA.(w)  ,  (7) 

J  u.'--o  J 

where  dA AJ)  is  a  random  increment  function  for  an  interval  du,. 

J 

Formation  of  the  u^  correlation  matrix  with  lag  i  gives: 
uj  i (  t)  =  <  Uj(t)  u^t+i)  > 

oo 

=  2  Re  /  exp(-iwi)  <  dA.(u»)dA-|(  <*>)  >  ,  (8) 

w=o  ^ 

where  *  denotes  complex  conjugate.  If  t  =  0,  U.,(  j  reduces  to  a 

J  * 

covariance  matrix.  The  covariance  matrix  formed  by  the  dA.(^)  gives 

J 


the  contribution  to  U  - , ( i)  by  the  cross  spectrum  i.  U)  in  an 

o  I  J  • 

infinitesimal  frequency  interval  d^: 

d'  =  \  <  dA j ( co)  dA* {  w)  >  .  (9) 

If  a  filtered  series  u.^(t)  is  defined 

J 

u.f(t)  =  Re[exp( iwt)  dA.(cj)]  du;  ,  (10) 

J  J 

then  the  filtered  correlation  matrix  U,  f(t)  is 

■J  • 

uj  1  f ( T )  =  Re[exp(-iu)i)$^((j)]  .  (11) 

If  PCA  is  applied  to  the  filtered  time  series,  it  would  involve  finding 
the  eigenstructure  of  the  symmetric  matrix  U.,f(o),  or  equivalently 

J  * 

Re[ci.(u))].  But  this  would  lose  the  phase  information  contained  in  the 
quadrature  components  of  the  i.ross  spectra,  which  is  the  very 
information  needed  to  detect  and  separate  traveling  waves  in  the  data 
set. 

Since  is  Hermitian,  it  has  a  complete  set  of  orthonormal 
eigenvectors  and  real  eigenvalues,  but  the  eigenvectors  are  now 
complex.  Wallace  and  Dickinson  (1972)  then  showed  how  to  find  this 
ei genstructure  which  was  then  applied  to  an  augmented  time  series, 
v(t),  involving  the  original  time  series  and  its  time  derivative,  so 
that  the  component  scores,  Z,  took  the  form 

z.(t)  =  Re  [  ei jVj(t)  , 

J 

where  e^j  is  the  appropriate  eigenvector. 


(12) 


d.  time  domain  CPCA 


While  FDPCA  would  allow  the  study  of  propagating  phenomena  in  a 
data  set,  it  can  be  cumbersome  in  certain  situations.  If  the  power  of 
a  component  is  spread  over  a  wide  frequency  band,  then  separate  maps 
for  each  spectral  estimate  must  be  studied  in  order  to  document  the 
propagating  phenomenon  (Horel,  1984).  Rasmusson  et  al.  (1981),  Barnett 
(1983),  and  Horel  (1984)  have  all  developed  and  documented  the  idea  of 
CPCA  in  the  time  domain  (hereafter  referred  to  simply  as  CPCA)  which 
can  be  easier  to  use  than  FDPCA.  CPCA  is  essentially  FDPCA  averaged 
over  all  frequency  bands.  CPCA  is  an  efficient  detector  of  propagating 
features,  especially  when  the  variance  is  spread  over  a  number  of 
frequencies  (Horel,  1984). 

Description  of  propagating  features  (waves)  in  a  data  set  requires 
both  co-spectrum  and  quadrature  information.  In  the  FDPCA  as  described 
by  Wallace  and  Dickinson  (1972),  this  information  was  contained  in  the 
cross-spectrum  matrix  and  extracted  and  contained  in  the 

ei genstructure  of  the  matrix  so  propagating  disturbances  could  be 
detected.  In  the  time  domain  this  information  is  not  so  readily 
available.  This  information  can  be  obtained  through  the  use  of  the 
Hilbert  transform.  The  following  theory  describing  the  use  of  the 
Hilbert  transform  is  based  on  the  work  of  Barnett  (1983). 

A  scalar  field  u(xj,t),  where  j  is  a  spatial  position  index  and  t 
is  time,  has  a  Fourier  representation  of  the  form: 


Uj( t)  =  .  a j (  • )  cos  t  +  bj( . )  sin  .t  , 


where  the  Fourier  coefficients  (a,b)  at  positive  frequencies  .  are 


defined  in  the  usual  manner.  The  variable,  u,  is  then  described  by  a 
complex  representation,  i.e., 

u.(t)  =  I  CjM  e_i  ,  (14) 

where  c-(uj)  =  a-{u>)  +  ib-(u)  and  i  =  ‘-1  .  Expanding  (14)  and 

vJ 

collecting  terms  yields 

Mt)  =  I  { [a j ( -■)  cos  ut  +  b-(u)  sin  4] 

u. 

+  i [ b j ( - )  cos  .-t  -  aj(-)  sin  .t]} 

=  Uj(t)  +  i  Uj(t)  .  (15) 

The  real  part  of  (15)  is  just  u-(t)  from  (13).  The  imaginary  part  of 

J 

(15),  u-(t),  represents  u-(t)  phase  shifted  ?/2  in  time  and  is  called 

J 

the  quadrature  function  or  Hilbert  transform  of  u-(t). 

J 

Wallace  and  Dickinson  (1972)  had  quadrature  information  available 
in  the  frequency  domain  from  (11)  unless  they  used  zero  lag,  in  which 
case  there  would  only  have  been  the  same  information  that  would  be 
available  from  a  real  valued  covariance  matrix  in  the  time  domain. 

This  zero  lag  covariance  matrix  problem  is  eliminated  in  the  time 
domain  by  going  into  the  frequency  domain  through  the  Fourier 
transform,  (13),  and  then  summing  over  frequencies  the  quadrature 
information  contained  there,  (15),  and  capturing  it  as  the  imaginary 
part  of  a  complex  number  when  returning  to  the  time  domain  via  the 
inverse  Fourier  transformation.  This  u  is  roughly  equivalent  in  the 
time  domain  to  the  quadrature  spectrum  in  the  frequency  domain.  This 
is  because  the  -/2  phase  shift  of  incremental  Fourier  component  is 


equivalent  to  providing  optimum  lag  information  on  that  component. 

Thus,  in  a  sense,  the  information  on  relative  phase  among  the  elemental 
sinusoids  has  been  summed  together.  It  is  also  obvious  that  the 
frequency  bandwidth  of  interest  (time  period)  can  be  limited  by 
applying  a  band-pass  filter  to  the  original  data  set,  or  by  limiting 
the  frequency  range  of  summation  in  (15). 

The  covariance  matrix  of  U.(t)  is 

SJ 

Cjk  =  <  Uj*(t)  Uk(t)  >  t  ’  (16) 

where  *  denotes  complex  conjugation  and  <  ...  indicates  a  time 
averaging  process.  Usually  would  be  a  zero  lag  covariance  matrix 
containing  no  phase  relationship  information.  But  in  this  case, 
has  been  constructed  from  U(t)  which  has  the  phase  information  built 
into  it.  Since  C  is  the  complex  "zero  lag"  covariance  matrix,  it  can 
be  shown  from  standard  spectral  analysis  theory  that 


where  ^  is  the  complex  cross-spectrum  and  the  Nyquist  frequency. 
So  represents  the  cross  spectral  matrix  averaged  over  a  specific 
frequency  band. 


C  is  Hermitian  by  construction  and  has  real  eigenvalues  (x  )  and 
complex  eigenvectors,  Bn(x),  where  x  represents  spatial  dependence.  U 
has  a  representation  which  optimally  accounts  for  the  variance  in 


u(x,t)  in  the  time  period  (frequency  band)  under  study,  i.e., 


where  the  complex  time-dependent  component  scores  are  given  by 

An(t)  =  I  U(x,t)  Bn(x)  .  (19) 

Four  measures  that  define  possible  moving  features  in  u(x,t)  can 
be  defined.  To  help  clarify  the  interpretation  of  these  measures  a 
simple  example  is  carried  along  with  each  definition  with  u(x,t)  =  a 
sin(kx-wot) ,  i.e.,  a  wave  of  constant  amplitude  and  frequency 
propagating  in  the  positive  x  direction.  The  measures  are: 

1)  spatial  phase  function  en(x):  this  function  shows  the 
relative  phase  of  fluctuations  among  the  various  spatial  locations 
where  u  is  defined.  This  measure  varies  continuously  between  0 
and  2 iT)  and  an  arbitrary  reference  value  must  be  selected  for  it. 
The  function  is  defined  as 

ep{x)  =  arctan  [Im  Bn(x)/Re  Bn(x) ]  .  (20) 

For  the  example  being  used,  the  data  field  is  composed  of  a  single 
sinusoid  with  wavelength  L  moving  along  the  x  axis.  There  is  only 
one  component,  and  e^x)  will  go  through  one  complete  cycle  (2tt) 
over  a  distance  L  measured  in  the  x  direction,  i.e.,  e^(x)  =  kx  = 
2ttx/L.  If  the  data  set  covered  a  spatial  area  of  8000  km  in  width 
in  the  x  direction  and  the  sinusoid  had  a  wavelength  of  4000  km, 
then  the  spatial  phase  function  should  progress  for  0  to  360 
degrees  twice  in  the  x  direction.  [Note  that  the  VAX/VMS  version 
4  FORTRAN  intrinsic  function  ATAN2D  was  used  in  the  computation  of 
the  spatial  and  temporal  phase  function,  and  it  has  a  range  from 
-180  <  result  <  180  degrees.] 


2)  spatial  amplitude  function  Sn ( x ) :  a  measure  of  the 
spatial  distribution  of  variability  associated  with  each 
component.  The  distribution  of  S  gives  a  measure  of  the  spatial 
homogeneity,  by  component,  in  the  magnitude  of  the  u-field. 

S„(x)  -  £Bn(x)  Bn*(x)]1/Z  .  (21) 

For  the  constant  amolitude  sinusoid  example  being  used,  equals 
a  constant  since  the  amplitude,  a,  is  independent  of  x.  The 
spatial  plot  of  this  function  would  show  equal  values  of  amplitude 
at  all  grid  points. 

3)  temporal  phase  function  <{>  (t):  this  measure  describes  the 
temporal  variation  of  phase  associated  with  u(x,t),  by  component, 
and  is  defined  as 

<t>n(t)  =  arctan  [Im  An(t)/Re  An(t)]  .  (22) 

The  example  sinusoid  of  fixed  frequency  and  wavenumber  will  have 
4> i ( t )  =  <*>  t  =  constant  x  t.  The  phase  will  progress  from  0  to 
360°  over  equal  time  intervals. 

4)  temporal  amplitude  function  Rn ( t ) :  this  measure  of 
temporal  variability  in  the  magnitude  of  the  modal  structure  of 
the  u-field  is  defined  as 

Rn(t)  -  [An(t)An*(t)]1/Z  .  (23) 

The  example  sinusoid  has  =  constant  over  time,  i.e.,  the  plot 
will  be  a  straight  line  at  a  constant  value. 

The  functions  described  in  (20)  to  (23)  above  constitute  a  general 
description  of  possible  moving  features  in  u.  Extraction  of  the 
components  of  the  covariance  matrix,  C-.  ,  allows  the  study  of  the 
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largest  spatially  coherent  "signals"  present  in  the  data  set.  This 


analysis  technique  allows  the  eas.,  investigation,  in  time  and  space,  of 
the  dynamic  behavior  of  energetic  events  that  might  well  be  obscured  in 
a  normal  spectral  analysis  approach  (Barnett,  1983). 

e.  methods 

The  data  set  used  in  this  study,  as  described  earlier,  consisted 
of  the  deviations  of  SLP  from  the  individual  grid  point  yearly  means. 
These  eight  sets  of  deviations  (180  observations  in  each  time  series  at 
280  grid  points)  were  filtered  to  concentrate  on  a  certain  frequency 
band  (time  period)  (described  in  section  2),  and  to  augment  the 
deviations  with  their  Hilbert  transform.  The  filter  used  was  developed 
from  the  theory  presented  by  Cizek  (1970).  The  filter  had  the  form: 

0  ,  K  =  1,2,. ..,135  and  173,  174 .  180 

H,  =  (24) 

K  +2  ,  K  =  136,  137,  ...,  172  . 

This  filter  was  applied  to  the  data  set  while  in  the  frequency  domain, 
(13),  and  produced  the  desired  band-passed  filtered  data  set  (2-10 
days)  with  its  Hilbert  transform  via  the  inverse  Fourier  transform. 

These  eight  yearly  filtered  data  sets  were  used  to  form  eight 
Hermitian  covariance  matrices  which  were  then  summed.  The  summed 
matrix  was  then  diagonalized  to  find  the  eigenvectors  and  eigenvalues 
using  the  IMSL  ( International  Mathematical  and  Statistical  Library) 
subroutine  EIGCH.  The  eigenvectors  (components)  were  ordered  by  their 
associated  eigenvalues,  and  then  the  component  scores  were  computed. 
From  the  components  and  component  scores,  the  four  functions  as 
described  in  section  3 . d . 1 ) -4 )  were  computed  for  the  four  components 


having  the  four  largest  eigenvalues.  These  four  components  were 
selected  for  study  after  preliminary  analysis  of  the  spatial  phase  and 
amplitude  function  plots  revealed  they  could  be  interpreted  as  synoptic 
events.  Components  5  through  280  were  not  interpretable,  and 
therefore,  for  this  study's  purposes,  considered  noise. 

The  spatial  amplitude  and  phase  plots  were  analyzed  using  a 
computer  graphics  program  which  put  the  spatial  phase  plot  "in  motion" 
while  overlaid  on  the  spatial  amplitude  plot.  A  color  scheme  was  used 
where  over  a  360-degree  phase  angle  change  the  color  representing  the 
phase  would  change  from  red,  to  green,  to  blue,  and  then  start  over 
again.  This  changing  of  colors,  as  the  computer  scans  the  phase  angle 
data  and  assigns  a  color  based  on  the  phase  angle  value,  appears  to  put 
the  plot  "in  motion."  This  greatly  enhances  the  ability  to  see  the 
motion  of  the  fluctuations  in  the  data,  as  represented  by  the  spatial 
phase  angle  changes.  The  values  of  the  spatial  amplitude  plot  were 
used  to  control  the  intensity  of  the  colors,  so  the  phase  changes  could 
be  detected  more  easily  in  the  areas  of  greater  variability.  Results 
are  reported  in  section  5. 


4.  INTERPRETATION  LIMITATIONS 


Before  considering  the  results  of  the  interpretation  of  the 
component  patterns  from  the  CPCA,  it  is  important  to  understand  some  of 
the  limitations  that  affect  the  interpretation.  Some  of  the 
limitations  can  be  overcome  through  rotation  of  the  initially  computed 
components,  and  these  techniques  have  been  demonstrated  by  Richman 
(1981)  and  Horel  (1981)  for  PCA.  Horel  (1984)  outlined  the  method  to 
be  used  for  CPCA,  but  it  appears  to  have  a  shortcoming,  which  is  also 
discussed  in  this  section. 

One  of  the  limitations  of  any  type  of  principal  component  analysis 
is  that  the  component  shapes  depend  on  the  geometry  of  the  domain, 
since  the  covariance  between  two  variables  can  be  considered  the  kernel 
of  an  integral  equation  eigenvalue  problem  (e.g.,  see  North  et  aK, 
1982).  The  integral  extends  over  the  finite  domain  of  interest.  So  if 
the  shape  or  size  of  the  domain  is  changed,  the  component  patterns 
should  be  expected  to  change. 

When  computing  the  components  of  the  data  set,  the  first  component 
is  positioned  in  hyperspace  to  account  for  the  largest  amount  of 
variance  over  all  the  iata  points,  with  each  succeeding  component 
orthogonal  to  it  and  to  the  other  components  to  follow.  Each 
succeeding  component  accounts  for  the  maximum  residual  variance 
remaining  in  the  input  data  matrix.  Prorlems  frequently  arise  when  the 
first  component  locates  between  two  or  more  distinct  clusters  of 
interrelated  variables.  In  this  study  these  variables  would  be 
stations  located  nearest  to  synoptic  cyclone  or  anticyclone  anomaly 
systems  (Kutzbach,  1967;  Richman,  1981).  Another  rotation  would  be 
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needed  to  accurately  locate  the  data  clusters  so  the  component  maps  can 
be  interpreted  as  an  individual  synoptic  map  type  (e.g.,  see  Richman, 
1981,  for  comparison  of  rotation  techniques).  Rotated  patterns 
generally  contain  somewhat  stronger  anomalies  covering  smaller  portions 
cf  the  domain  of  interest.  This  "regionality"  of  fluctuations  is  very 
desirable. 

Due  to  the  orthogonality  constraint  placed  upon  the  components  of 
the  data  matrix,  they  may  take  on  predictable  geometric  relationships 
to  one  another.  If  the  first  component  is  of  one  sign  over  the  domain 
(applies  to  PCA),  the  second  component  pattern  may  be  bipolar  and  have 
a  zero  line  passing  through  the  highest  value  of  the  first  component 
pattern,  and  so  on  for  the  other  components  (Horel ,  1981).  In  the  case 
of  CPCA,  all  component  variability  patterns  are  positively  valued,  but 
low  variability  areas  on  one  map  may  match  the  high  variability  areas 
on  another. 

Linear  transformation  (rotation)  of  the  components  is  performed  to 
overcome  the  problems  mentioned  above.  A  variety  of  procedures  exist 
to  accomplish  the  rotation,  e.g.,  varimax,  quartimax,  promax,  oblimax 
(Richman,  1981;  Horel,  1981;  Green,  1978),  all  of  which  should  identify 
regional  relationships  from  the  initial  component  solution.  In  the 
case  of  conventional  PCA,  using  real  data  sets,  these  rotation 
procedures  present  no  problems. 

In  the  PCA  of  a  data  set  the  components  can  be  thought  of  as 
orthogonal  lines  that  bisect  the  axes'  origin.  These  lines  have  both  a 
positive  and  negative  valued  segment.  Either  segment  can  be  used  as  a 
measure  of  the  variance  accounted  for  by  the  component  of  interest. 
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The  negative  or  positive  coefficient  is  associated  with  the  phase 
(direction  of  the  line)  of  the  component.  The  phase  is  only  known  to 
within  an  additive  constant  (either  0  or  i),  so  the  sign  of  the 
component  is  ambiguous. 

In  the  CPCA  approach  the  phase  angle  is  known  only  to  within  an 
additive  constant  that  ranges  from  0  to  2u.  This  indeterminancy  is  a 
problem  when  linearly  transforming  (rotating)  the  complex  components. 


As  outlined  by  Horel  (1984),  using  the  varimax  rotation  technique,  an 
element,  e-  ,  of  the  nth  complex  eigenvector  of  a  complex  correlation 

J  1  * 

matrix,  can  be  written  as 


ejn  *  <  >t  ’  V’6j" 


(25) 


where  e-  is  the  complex  correlation  between  the  jth  time  series  and 

J  *  * 

the  nth  component,  Sjn  is  the  magnitude  of  the  correlation,  and  is 
the  phase.  An  estimate  UUt)  of  the  observation  matrix  U  - ( t )  can  be 

J  J 

reconstructed  with  the  first  N'  complex  components  where  N'  <  N,  where 
N  is  the  total  number  of  complex  components  as 


N' 


Uj(t)  '  Jj  ejnAn(t)  • 


(26) 


The  estimate  UUt)  can  be  described  in  terms  of  the  transformed 

J 

(rotated)  solution  as 
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uj(t>  ■  j,  iF.<‘»  • 


(27) 


where  F  (t)  is  an  element  of  the  nth  rotated  complex  correlation  (RCi  C) 
m 

and  gjm  is  the  complex  correlation  between  the  nth  RCPC  and  the  jth 
observation  of  IK(t) . 

J 

In  PCA  the  total  amount  of  variance  explained  by  the  first 
component  is  maximized 


U  *\ 

y  e? 
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(28) 


where  e^  is  the  (real)  correlation  between  the  jth  time  series  and  the 
first  component.  Likewise,  CPCA  maximized 


V/ 
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where  e^  is  a  complex  correlation  and  e^  is  its  complex  conjugate. 
PCA  and  CPCA  maximize  the  sum  of  the  squares  of  the  amplitudes  of  the 
correlation  coefficients.  The  varimax  rotation  technique  maximizes  a 
second  moment  statistic  (the  variance),  rather  than  a  first  moment 
statistic  (the  sum)  of  the  squared  correlation  coefficient.  This 
causes  strong  regional  relationships  in  large  data  sets  to  be 
emphasized. 

In  order  to  achieve  this  variance  maximization,  the  complex 


1 


loadings  gjm  are  expressed  in  terms  of  the  complex  eigenvectors  e j n : 


yjm  l  mn  jn 
n=l 

where  the  coefficient  t _  of  the  transformation  (rotation)  matrix  T  are 

real.  The  maximization  of  the  variance  is  dependent  on  e-n,  which  can 
be  seen  from  (25),  to  depend  on  a  phase  angle  known  only  to  within  an 
arbitrary  additive  constant.  Unlike  the  PCA  case  where  the  phase 
angles  are  0  or  it,  the  CPCA  has  as  possible  solutions  a  whole  set  of 
vectors  ranging  from  0  to  2 n,  each  as  valid  a  solution  as  the  other. 
Thus,  it  appears  that  a  maximization  of  the  variance  for  a  whole  set  of 
ejn's  (as  a  function  of  the  phase  angle)  is  needed  and  then 
determination  of  which  one  gives  the  "best"  maximization  of  the 
variance.  If  two  or  more  maximizations  are  nearly  equal,  then  a 
selection  of  the  best  solution  must  be  based  on  the  physical 
interpretation  of  the  solutions--which  one  gives  the  best  physical 
representation  of  the  data. 

Therefore,  the  initial  components  found  in  this  study  were  not 
rotated,  but  rather  interpreted  as  they  were.  More  work  needs  to  be 
done  in  the  area  of  rotation  of  complex  components,  so  that  the 
rotation  will  not  be  dependent  on  an  arbitrary  phase  angle.  Then  the 
rotation  procedure  will  be  of  a  manageable  nature,  where  it  doesn't 
require  the  computation  of  a  large  number  of  variance  maximizations. 

To  see  if  the  four  unrotated  components  studied  were  suffering 


from  the  limitations  outlined  above,  a  squared  coherence  check  was 
performed  between  the  components.  The  components  are  orthogonal  at 
zero  lag,  but  may  have  coherence  at  any  frequency.  The  squared 


coherence  check  reveals  if  there  is  coherence  present.  Table  1  shows 
the  squared  coherence,  by  frequency,  for  the  pairs  of  components.  A 
value  of  .122  equals  the  95  percent  point  of  the  null  distribution  of 
squared  coherence.  The  99  percent  point  equals  .181. 

The  spectra  of  the  four  components,  figure  1,  show  the  important 
periods  to  be  6.4  to  5.6  and  5.0  to  4.5  days  for  component  one;  9.0  to 
8.2  days  for  component  two;  9.0  to  6.4  days  for  component  three;  and 
9.0  and  6.9  days  for  component  four.  An  examination  of  table  1  shows 
that  all  pairs  of  components,  except  2  and  4,  have  significant 
coherence  at  the  periods  of  importance.  Thus,  we  can  assume  that  the 
individual  components  are  not  defining  distinct  clusters  of  variables 
(stations)  that  are  associated  with  synoptic  anomaly  systems.  They 
share  variance  (except  between  2  and  4)  from  the  same  regional  anomaly 
areas  over  the  domain,  instead  of  each  one  defining  one  distinct 
pattern  of  variance. 

The  spectra  studied  to  find  the  periods  listed  above  were 
constructed  by  summing  the  eight  yearly  periodograms  for  each 
component,  and  then  averaging  the  summed  periodogram  values  across 
three  frequencies  of  each  component's  periodogram.  This  resulted  in  48 
degrees  of  freedom.  The  95  and  99  percent  points  of  the  null 
distribution  of  squared  coherence  were  computed  via  the  equations 


Table  1.  Coherence  check  between  pairs  of  components  being 

interpreted.  Results  are  significant  at  the  95.  level  of 

o 

confidence  for  (coherence)  =  0.122 


p 

(coherence) 


Component  pairs:  1  vs  2  1  vs  3  1  vs  4  2  vs  3  2  vs  4  3  vs  4 


31 


5.  RESULTS 


a.  mean  SLP  and  standard  deviation 

Figure  2a  of  unfiltered  average  SLP  is  very  similar  to  the  figure 
la  of  the  same  variable  by  Blackmon  et  al_.  (1977)  which  covers  the 
years  1963  to  1972  for  15  November  to  14  March.  Both  figures  mentioned 
above  also  compare  favorably  to  the  mean  SLP  maps  in  the  Normal  Weather 
Charts  for  the  Northern  Hemisphere. 

Figure  2b  of  the  standard  deviation  of  band-passed  SLP  is  very 
similar  to  the  figure  2b  of  the  same  variable  by  Blackmon  et  aK  (1977) 
for  the  same  time  period  as  mentioned  above.  A  comparison  of  the 
standard  deviation  plots  with  the  variance  of  the  band-passed  500  mb 
heights  (e.g.,  see  figure  5a  and  c  of  Blackmon,  1976),  standard 
deviation  of  the  band-passed  300  mb  height  ( e.g .,  see  figure  3b  of 
Blackmon  et  aK ,  1977),  and  the  average  wintertime  500  mb  vorticity 
(e.g.,  see  figure  11  of  Blackmon,  1976)  reveals  some  interesting 
meteorological  connections  between  the  different  levels.  All  the 
graphs  have  maximum  values  in  the  area  of  Nova  Scotia  and  Newfoundland, 
with  a  relative  secondary  maximum  in  the  area  of  the  Colorado 
cyclogenesis  region  in  the  average  vorticity  field  of  Blackmon's  figure 
11  (1976).  Synoptic  experience  and  past  studies  (Blackmon  et  al . , 
1984a;  Klein,  1967)  have  indicated  that  the  "average"  cyclone  tracks 
are  determined  to  a  large  extent  by  the  upper  level  general  circulation 
(particularly  700  mb),  and  should  produce  a  maximum  in  the  pressure 
deviations  and  maximum  vorticity  at  all  levels  between  40°  and  50°N 
over  eastern  North  America  as  the  general  circulation  would  tend  to 
guide  cyclones  to  this  area.  These  cyclones  are  usually  mature,  deep. 
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and  large  area  covering  synoptic  scale  disturbances.  Therefore,  the  | 

effects  of  these  cyclones  should  be  seen  from  the  surface  to  the  upper 
levels  of  the  atmosphere  (Hartmann,  1974),  while  the  secondary  maximum 

i 

in  the  500  mb  average  vorticity  field  (e.g.,  see  figure  11  of  Blackmon,  j 

1976)  is  probably  induced  by  the  combination  of  the  lee-side  cyclonic 
flow  and  the  maximum  cyclogenesis  cases  per  unit  area  occurring  in  that 
area  (Whittaker  and  Horn,  1981).  j 

When  looking  only  at  the  band-passed  SLP  standard  deviations 
(figure  2b)  the  enhanced  SLP  fluctuations  can  be  seen  along  the  east  ; 

side  of  the  major  mountain  ranges  caused  by  cyclonic  flow  on  the  j 

lee-side  of  the  ranges.  Note  how  the  band-passed  SLP  standard  ] 
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deviation  contours  "pinch  in"  on  the  main  axis  of  deviations,  from  the  ! 

north  and  south,  from  approximately  100°  to  85°W.  Storms  moving  out  of  j 

the  lee-side  cyclogenesis  area  are  usually  weak  and  small  in  terms  of 
area  covered,  which  will  tend  to  cause  smaller  deviations  in  the  SLP 

t 

field,  than  when  they  initially  had  the  lee-side  circulation  support.  | 

I 

As  they  continue  to  cross  the  central  portion  of  the  United  States  or 

t 

southern  Canada  the  cyclones  tap  moisture  and  heat  sources  (Gulf  of  j 

i 

Mexico  or  Great  Lakes)  and  begin  to  intensify,  covering  larger  areas  of  I 

the  earth,  and  therefore  widen  and  strengthen  the  standard  deviation  i 

1 

field  of  the  band-passed  SLP.  Results  of  Hartmann  (1974),  Blackmon  et  j 

al .  (1977,  1979),  and  Blackmon  and  White  (1982)  support  the  correlation  j 

between  the  standard  deviation  field  in  SLP  and  upper  level  atmospheric 

parameters  when  moving  from  the  highly  baroclinic  disturbance  areas  of 

lee-side  cyclogenesis  to  the  less  baroclinic,  more  mature  cyclones  off  j 

( 

the  east  coast  of  northern  North  America.  < 


As  noted  by  Blackmon  et  a^.  (1977),  a  comparison  of  band-passed 


SLP  standard  deviations  (figure  2b)  with  the  statistics  derived  by 
Petterssen  (1956)  of  the  percentage  frequency  of  cyclone  centers  per 
unit  area  during  the  winter  season  provides  some  interesting 
information.  Petterssen's  figures  13.6.1  and  13.6.2  show  maxima  of 
cyclone  centers  over  the  Great  Lakes  and  other  warm  bodies  of  water, 
whereas  this  study's  results  and  those  of  Blackmon  et  aK  (1977) 
indicate  only  a  little  enhancement  of  pressure  variability  in  these 
areas.  This  suggests  that  many  of  Petterssen's  centers  over  these 
areas  are  weak  and/or  quasi-stationary  (Blackmon  et  aK ,  1977).  This 
would  be  analogous  to  the  "pinching"  seen  in  the  SLP  deviation  field  as 
disturbances  move  out  of  the  lee-side  cyclogenesis  areas  as  explained 
above. 

Blackmon  et  aK  (1977)  further  felt  that  their  variance  plot 
indicated  a  more  widespread  region  of  "lee-side  cyclogenesis" 
downstream  from  the  Rockies.  It  appears  that  the  band-passed  variance 
plots  indicate  the  movement  of  cyclones  out  of  the  lee-side  area  as 
weak  cyclones,  causing  "pinching"  in  the  SLP  deviation  field  before 
they  begin  to  later  intensify  and  expand  in  area  as  they  head  toward 
the  east  coast  of  North  America.  As  stated  by  Blackmon  et  aK  (1977), 
the  work  of  Reitan  (1974)  closely  resembles  the  band-passed  SLP 
standard  deviation  field,  which  it  should  since  the  cyclones  tend  to 
follow  the  longwave  circulation  and  trace  out  a  pattern  of  increasing 
pressure  variance  as  they  intensify  and  grow  in  areal  size.  Reitan's 
(1974)  figure  12  further  confirms  that  the  cyclone  tracks  to  the 
lee-side  of  the  Rockies  all  converge  in  the  Nova  Scotia-Newfoundland 
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area,  as  would  be  caused  by  the  mean  upper  level  flow  (e.g.,  see  figure 
2  of  Blackmon  et  ah  ,  1984b). 


The  results  of  Blackmon  (1976),  Blackmon  et  al .  (1977,  1979, 
1984a),  and  others  mentioned  in  this  section  agree  with  the  average  and 
standard  deviation  of  SLP  plots  (figures  2a  and  b)  of  this  study.  It 
doesn't  appear  the  data  set  being  considered  in  this  study  is  abnormal. 
The  results  of  the  CPC A  should  give  a  good  idea  of  the  normal  areas  of 
preferred  synoptic  scale  disturbances  and  their  movement. 

b.  component  number  one 

The  first  component  of  the  CPCA  accounted  for  approximately  26.6 
percent  of  the  variance  in  the  band-passed  filtered  data  set.  The 
period  of  the  most  important  wave  structures  contributing  to  the 
spectrum  of  this  component  were  estimated  to  be  4.5  to  5.0  days  and  5.6 
to  6.4  days. 

1)  spatial  properties 

Figure  3a  is  the  spatial  amplitude  plot  of  the  first  component. 
There  is  one  main  area  of  variability  centered  over  northern  Nova 
Scotia  with  a  ridge  of  variability  extending  to  the  NE  along  50°N,  and 
also  extending  to  the  Sin/  to  a  secondary  maximum  over  northern  Oklahoma. 
There  are  two  areas  of  lower  variability  evident,  the  first  being  the 
area  over  the  Atlantic  Ocean  near  21°N,  50°W,  which  should  be 
associated  with  the  subtropical  anticyclone  in  this  area.  There  is  an 
extension  of  this  area  along  the  Gulf  Coast  and  up  the  west  side  of  the 
Appalachian  Mountains.  The  second  area  of  lower  variability  is  located 
over  the  western  Hudson  Bay  area  extending  down  to  the  western  Great 
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Lakes.  This  area  appears  to  be  induced  between  the  higher  variability 
areas  of  the  lee-side  of  the  Rockies  and  Nova  Scotia. 

The  spatial  amplitude  plot  closely  resembles  the  average 
wintertime  500  mb  vorticity  pattern  (e.g.,  see  figure  11  of  Blackmon, 
1976).  The  area  of  maximum  values  in  both  plots  over  Nova  Scotia 
agrees  well  with  the  idea  that  cyclones  that  reach  this  area  are 
mature,  deep  disturbances  affecting  much  of  the  atmosphere  vertically 
(Blackmon  et  aK,  1979;  Hartmann,  1974).  A  look  at  the  plots  of  500  mb 
geopotential  height  variance  (e.g.,  see  figure  5c  of  Blackmon,  1976) 
and  300  mb  height  standard  deviations  (e.g.,  see  figure  3b  of  Blackmon 
et  aK,  1977)  of  band-passed  data  also  observational  ly  supports  the 
idea  of  mature  cyclone  structures  near  Nova  Scotia. 

The  area  of  maximum  variability  in  figure  3a  is  also  located  to 
the  north  side  of  the  winter  average  500  mb  zonal  wind  speed  maximum 
(e.g.,  see  figure  6a  of  Blackmon  et  aj_. ,  1977).  This  location  is  in 
agreement  with  Palmen  and  Newton  (1969),  who  found  that  as  cyclones 
mature  the  sea  level  circulation  center  deviates  increasingly  to  the 
left  of  the  jet  stream. 

White  (1980)  plotted  the  ratio  of  band-passed  wintertime  500  mb 
geopotential  height  variance  contribution  of  long  waves  to 
synoptic  scale  waves.  White's  figure  2  shows  a  distinct  minimum  in  the 
ratio  over  the  Nova  Scotia-Newfoundland  area,  where  component  one  has 
maximum  variability  due  to  synoptic  scale  activity.  This  ratio  minimum 
shows  that  the  long  waves  are  contributing  very  little  to  the  variance 
in  the  500  mb  geopotential  height,  further  supporting  the  importance  of 
synoptic  scale  activity  in  the  Nova  Scotia-Newfoundland  area. 


Figure  3b  is  the  spatial  phase  plot  for  component  one.  The 
relative  phase  angles  increase  from  west  to  east,  indicating  eastward 
movement  of  the  fluctuations  in  the  SLP  field.  The  more  widely  spaced 
isolines  of  phase  angles  (in  units  of  degrees)  in  the  lee-side  of  the 
Rockies  indicates  the  faster  movement  of  these  young,  weak  disturbances 
as  they  move  in  a  southern  direction  before  turning  eastward.  After 
turning  eastward  the  phase  angle  isolines  in  the  northern  portion  of 
the  United  States  indicate  a  relative  decrease  in  the  area  covered 
between  changes  in  isoline  values.  This  appears  to  indicate  a  relative 
speed  decrease  in  the  movement  of  the  disturbances  in  the  northern 
United  States.  Meanwhile,  in  the  southern  portion  of  the  United  States 
the  phase  angle  isolines  appear  to  indicate  the  relative  speed  of  the 
disturbances  has  not  slowed  down  as  much  as  they  have  in  the  northern 
United  States,  until  they  start  being  influenced  (blocked)  by  the 
Appalachian  Mountains.  The  disturbances  in  the  northern  United  States 
appear  to  speed  up  once  they  reach  the  Great  Lakes  area,  and  all 
disturbances  appear  to  increase  in  speed  once  they  come  under  the 
influence  of  the  warm  Atlantic  Ocean  waters,  especially  in  the  area  of 
the  Gulf  Stream.  One  possible  explanation  for  this  speed  increase 
would  be  that  the  disturbances  are  receiving  energy  from  the  warm 
waters  of  the  Great  Lakes,  Gulf  of  Mexico,  or  Atlantic  Ocean.  More 
work  needs  to  be  done  to  support  this  explanation.  The  question  to  be 
answered  is,  do  disturbances  in  fact  speed  up,  on  a  consistent  basis, 
when  energy  is  input?  It  also  appears  from  the  plot  that  once  the 
disturbances  move  east  of  60°W,  they  move  at  a  fairly  regular  speed. 
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These  disturbances  should  be  mature,  fully  developed  systems,  which 
have  reached  a  fairly  steady  state  of  energy  input. 

The  structure  of  the  phase  angle  isolines  in  figure  3b  correspond 
closely  to  the  structure  of  the  patterns  of  correlation  values  of  500 
mb  height  deviations  which  exhibit  a  west  to  east  oriented  wave-like 
structure,  with  a  zonal  wavelength  of  about  4000  km,  as  reported  by 
Blackmon  et  aj_.  (1984a).  The  results  of  Blackmon  et  (1984a) 
indicate  that  these  short-time  scale  (2.5  to  6.0  days)  fluctuations  are 
aligned  slightly  poleward  and  downstream  of  the  climatological  mean  jet 
stream  speed  maximum  and  with  axes  oriented  along  the  climatological 
mean  flow.  The  fluctuations  which  produce  these  correlation  patterns 
also  have  a  predominant  zonal  scale  of  wavenumber  6-7,  and  are  clearly 
related  to  baroclinic  eddies  (synoptic  scale  disturbances)  as  discussed 
by  Blackmon  M976)  and  Blackmon  et  aK  (1977). 

It  appears  the  700-500  mb  mean  wind  flow  pattern  is  creating  the 
phase  angle  pattern  shown  in  figure  3b.  An  examination  of  the  mean 
wintertime  700  and  500  mb  flow,  and  the  results  of  Klein  (1967)  and 
Blackmon  et  al_.  (1984b)  all  support  the  idea  that  the  "average" 
variance  pattern  in  SLP  is  traced  out  by  the  upper  level  mean  flow 
pattern.  Further  support  is  given  by  the  work  of  Bowie  and  Weightman 
(1914),  Reitan  (1974),  and  others  showing  the  convergence  of  preferred 
storm  tracks  into  the  area  near  Nova  Scotia  after  having  followed  paths 
across  North  America  that  approximate  the  mean  upper  level  flow 
pattern. 


2)  temporal  properties 

Figure  3c  is  the  temporal  amplitude  plot,  by  year,  for  component 
one.  The  plot  shows  that  component  one  was  weakest  in  1981  and 
strongest  in  1973.  This  can  be  seen  by  looking  at  the  squared  temporal 
amplitude  (variance)  in  table  2  also. 

The  monthly  weather  surmiaries  in  Weatherwise  and  the 
Climatological  National  Summary  suggest  that  1973  was  the  most  active 
year  (of  the  eight  studied)  of  cyclone  activity,  with  the  700  mb  jet 
displaced  south  to  the  gulf  coast.  On  the  other  hand,  1981  was  a  year 
with  a  blocking  ridge  along  the  west  coast  and  a  trough  along  the  east 
coast  causing  a  large  number  of  storms  to  form  and  intensify  as  they 
moved  off  the  coast  of  Nova  Scotia.  Therefore,  they  did  not  contribute 
a  great  deal  to  the  "overall"  pattern  of  component  one. 

Figure  3d  is  the  temporal  phase  plot  of  component  number  one  which 
shows  that  it  is  dominated  by  disturbances  whose  phase  changes  with  a 
fairly  quick,  linear  progression.  Only  a  few  stand  out  as  being  slower 
progressing  features.  But  these  did  not  contribute  significantly  to 
the  variance  of  the  component  as  can  be  seen  by  comparing  the  temporal 
amplitude  plot  to  the  temporal  phase  plot  at  the  time  of  occurrence  of 
the  more  slowly  progressing  features. 

3)  comments 

Component  one  seems  to  represent  the  pattern  produced  in  the 
pressure  variance  field  by  disturbances  moving  roughly  with  the  mean 
mid-tropospheric  flow.  Care  must  be  taken  that  the  variance  is  not 
described  only  in  terms  of  cyclone  activity.  Petterssen  (195C) 
describes  the  "cyclone  family"  as  consisting  of  both  a  train  of 


Figure  3d.  The  temporal  phase  (degrees)  of  component  one  of  the 
band-passed  filtered  data,  by  year 
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Table  2.  Squared  temporal  amplitude  (variance),  by  year,  for  each 
component  in  units  of  squared  millibars. 


2 

variance  (mb) 


Year: 

73 

74 

75 

76 

79 

80 

81 

82 

Comp.  1 

7155 

5279 

5146 

6393 

5458 

5289 

4286 

6012 

Comp.  2 

3317 

3812 

3138 

4490 

2916 

2730 

2863 

3491 

Comp.  3 

3536 

2387 

2516 

3775 

2405 

2515 

2582 

3434 

Comp.  4 

1622 

2203 

1434 

1649 

1711 

1578 

709 

1479 

cyclones  and  anticyclones  (e.g.,  see  figure  12.5.1  and  12.5.2  of 
Petterssen,  1956).  The  anticyclones,  as  noted  by  Petterssen  (1956), 
appear  quite  often  to  be  sluggish,  passive  systems  that  fill  the  space 
between  active  cyclone  systems.  This  seems  to  fit  in  well  with  the 
occasional  more  slowly  changing  phase  angle  of  some  events  in  the 
temporal  phase  plot.  Klein  (1958)  found  that  the  average  life  span  for 
all  pressure  centers  was  approximately  five  days,  with  anticyclones 
usually  lasting  one  day  longer  than  cyclones. 

c.  component  number  two 

The  second  component  accounts  for  approximately  15.8  percent  of 
the  variance  in  the  band-passed  data  set.  The  period  of  the  most 
important  wave  structure  contributing  to  the  spectrum  of  this  component 
was  estimated  to  be  8.2  to  9.0  days. 

1)  spatial  properties 

Figure  4a  is  the  spatial  amplitude  plot  for  component  two.  There 
is  one  maximum  area  of  variability  over  the  Dakotas  and  Minnesota  with 
extensions  of  this  area  along  ridges  of  variability  southward  to  the 
Texas  gulf  coast,  northward  into  Alberta,  and  ENE  through  Canada  along 
and  just  north  of  50°N  latitude.  There  appears  to  be  a  relative 
secondary  maximum  in  the  variability  off  the  Atlantic  coast  centered 
near  35°N,  64°W.  The  two  areas  of  relative  minimum  variability  are 
located  over  Florida,  and  near  43°N,  53°W. 

The  spatial  amplitude  pattern  maximum  over  the  Dakotas  and 


Minnesota  matches  closely  the  area  of  maximum  frequency  of  anticyclones 
centers  found  by  Miller  (1932),  Zishka  and  Smith  (1980),  and  Klein 
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(1957).  The  extensions  of  this  area  to  the  south  and  ENE  follow 
closely  the  preferred  tracks  of  anticyclones  as  reported  by  Zishka  and 
Smith  (1980),  Wasko  (1954),  and  Klein  (1957). 

A  look  at  the  actual  cyclone  and  anticyclone  tracks,  for  the  eight 
years  studied,  shows  the  area  along  the  United  States-Canadian  border 
was  a  very  frequent  path  used  by  these  disturbances.  Further,  there 
were  a  number  of  anticyclones  that  moved  southward  into  the  central 
portion  of  the  United  States. 

An  examination  of  figure  9b  of  Blackmon  et  aK  (1977)  of 
band-passed  850  mb  temperature  standard  deviation  reveals  that  this 
plot  and  figure  4a  of  this  study  are  nearly  identical  in  shape.  The 
area  of  maximum  variability  in  the  SLP  to  the  east  of  the  Rocky 
Mountains  and  the  area  of  maximum  deviation  in  the  850  mb  temperature 
plot  of  Blackmon  et  al_.  (1977)  could  possibly  be  linked  to  the  outbreak 
of  shallow,  modified,  polar  air  into  the  central  United  States,  just 
east  of  the  Rockies,  during  the  winter  season.  These  cold  air 
outbreaks  can,  and  do,  reach  as  far  south  as  Texas  and  Mexico  (Blackmon 
et  aK ,  1979).  The  graphs  of  Blackmon  et  al_.  (1979)  of  the  correlation 
between  the  1000  mb  and  500  mb  height  (their  figure  3b)  and  between 
1000  mb  height  and  the  1000-500  mb  thickness  (their  figure  12)  supports 
the  idea  that  the  outbreaks  are  relatively  shallow  and  do  not  extend  up 
to  the  500  mb  level  in  the  atmosphere.  The  increase  in  pressure  at  the 
surface  is  compensated  for  by  the  decrease  in  temperature  in  the  lower 
atmosphere,  so  the  effects  of  the  polar  air  are  not  seen  at  the  upper 
levels  of  the  atmosphere.  A  look  at  surface  temperature  anomalies  for 
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the  period  of  study  shows  departures  from  normal  as  much  as  12  to  18 
degrees  Celsius  below  normal  over  the  west-central  United  States. 

The  area  of  low  variability  over  Florida  is  probably  present  due 
to  the  subtropical  high  being  in  control  of  the  area  most  of  the  time. 
The  other  low  variability  area  is  probably  induced  by  the  CPCA 
technique.  As  mentioned  in  section  4,  since  component  one  had  a  large 
variability  center  over  Nova  Scotia,  the  low  variability  at  43°N,  5  3°W 
in  component  two  could  be  induced  by  component  one's  pattern. 

Figure  4b  of  the  spatial  phase  plot  for  component  tw<  sh  ^  •*.» 

movement  of  the  fluctuations  has  two  main  paths.  Both  pat*-,  a;.-;  ;>  *. 
start  in  Canada  around  58°N,  116°W  and  travel  SSE  toward  the  .'  Te¬ 
states.  One  path  continues  into  the  United  States,  heading  toward 
Texas  and  the  gulf  coast,  while  the  other  breaks  off  and  heads  east, 
along  the  50°N  latitude  circle.  The  path  into  the  United  States  and 
toward  Texas  matches  well  with  the  idea  of  cold  air  outbreaks.  The 
other  path,  along  the  50°N  latitude  circle,  matches  with  the  idea  of  a 
major  anticyclone/cyclone  path  through  this  area  as  shown  by  Klein 
(1958)  and  Wasko  (1954). 

2)  temporal  properties 

The  temporal  amplitude  plot,  figure  4c,  for  component  two  shows 
that  1976  was  the  year  in  which  the  contribution  of  the  second 
component  to  the  total  variance  was  a  maximum,  while  1980  contributed 
the  least  (see  table  2  also).  In  1976  the  mean  700  mb  trough  was 
shifted  westward  and  centered  over  the  Great  Lakes.  This  shifted 
orientation  would  tend  to  direct  disturbances  southward  on  the  west 
side  of  the  trough  through  the  Great  Plains  area  more  than  usual. 
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While  in  1980,  the  west-central  United  States  experienced  a  warmer  than 
usual  winter.  There  were  few  cold  air  outbreaks  moving  southward  out 
of  Canada. 

The  temporal  phase  plot  (figure  4d)  reveals  that  the  phase  is 
changing  more  slowly  with  time,  for  the  most  part,  than  it  did  for 
component  one  (figure  3d).  This  is  in  keeping  with  the  idea  of  more 
slowly  moving  anticyclones  dominating  this  component  with  some  slower 
moving  cyclones  probably  mixed  in. 

3)  comments 

This  component  seems  to  be  associated  with  cyclones  which  are 
taking  a  northern  track  due  to  possibly  a  strong  anticyclone  over  the 
eastern  United  States  and  anticyclones  that  track  along  the  United 
States-Canadian  border.  It  is  also  linked  to  the  cold  air  outbreak 
anticyclones  that  track  southward  just  to  the  east  of  the  Rockies. 

These  outbreaks  seem  to  be  enhanced  by  a  trough  in  the  upper  level 
circulation  over  the  Great  Lakes. 

The  coherence  check  (table  1)  shows  that  this  component  is 
correlated  with  component  one,  especially  in  the  8.2  to  9.0  day  period. 
This  seems  to  indicate  the  cold  air  outbreaks  may  add  to  the 
variability  of  both  components,  just  to  the  east  of  the  Rockies,  due  to 
the  time  length  of  the  significant  period  which  would  indicate  slower 
moving  anticyclones. 


component  number  three 


The  third  component  accounts  approximately  for  13.7  percent  of  the 


variance  in  the  band-passed  data  set.  The  periods  of  the  most 


temporal  phase  (degrees 
!-passed  filtered  data, 
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important  wave  structures  contributing  to  the  spectrum  of  this 
component  were  estimated  to  be  6.4  to  9.0  days. 

1)  spatial  properties 

Figure  5a  is  the  spatial  amplitude  plot  for  the  third  component. 
There  are  two  areas  of  maximum  variability  separated  by  a  line  of  low 
variability  through  the  central  plains  states. 

The  area  of  variability  over  the  Rocky  Mountains  extends  out  of 
Alberta,  Canada,  down  to  the  Gulf  of  Mexico,  but  has  its  center  in  the 
central  Rockies  over  northern  Utah  and  southwest  Wyoming.  This  area 
corresponds  very  closely  to  a  maximum  area  of  anticyclone  events  and 
anticyclogenesis  as  depicted  by  Petterssen  (1956,  figures  13.7.1  and 
13.7.2)  and  Zishka  and  SmiJh  (1980).  The  pattern  is  also  located  in 
one  of  the  areas  of  high  alternation  between  cyclones  and  anticyclones 
in  the  winter  season  (e.g.,  see  figure  13.8.2  of  Petterssen,  1956). 
Petterssen  (1956)  speaks  of  polar  outbreak  highs  moving  along  the  west 
coast.  This  component  seems  to  represent  these  outbreaks  somewhat, 
while  component  two  seems  to  represent  the  polar-outbreak  highs  over 
the  central  plains  states. 

The  second  area  of  variability  centered  over  the  Coast  of  Labrador 
has  two  southward  curving  extensions:  one  extends  to  the  SW  to  the 
Great  Lakes,  and  the  other  one  extends  to  the  SE  to  49°N,  47°W,  and 
from  there  it  curves  back  to  the  SW  through  40°N,  50°W.  This  area 
seems  to  fit  the  cyclone  tracks  depicted  by  Reitan  (1974)  and  Zishka 
and  Smith  (1980)  for  the  area  extending  back  to  the  Great  Lakes  from 
the  center  over  the  Coast  of  Labrador.  But  the  area  extending  off  the 
Atlantic  Coast  appears  to  have  too  much  curvature  back  to  the  NW  after 
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passing  the  tip  of  Newfoundland,  to  match  very  well  with  published 
cyclone  tracks. 

The  extension  of  the  area  of  variability  off  the  Atlantic  Coast 
seems  to  be  artificially  forced  to  curve  due  to  the  induced  area  of  low 
variability  near  42°N,  68°W.  This  low  variability  area  appears  to  be 
present  due  to  the  extremely  high  variability  area  near  this  location 
in  component  one.  Since  component  one  accounted  for  a  large  part  of 
the  variability  in  this  area,  component  three  has  low  variability  there 
because  component  three  accounts  for  the  variability  left  in  the  data 
field  after  the  variability  due  to  components  one  and  two  has  been 
taken  out.  There  is  also  the  possibility  that  the  shape  of  the 
variability  area  is  caused  by  the  domain  problem  as  noted  in  section 
four.  Climatological  storm  track  data  show  that  one  of  the  primary 
tracks  for  storms  leaving  the  Great  Lakes  area  is  to  the  NE  through 
Quebec  and  the  St.  Lawrence  River  area.  At  the  same  time,  storms 
moving  around  the  eastern  side  of  Nova  Scotia  will  join  the  Great  Lakes 
storms  in  the  area  off  the  Labrador  Coast  and  then  all  of  them  move  off 
to  the  east  of  Greenland.  The  boundary  of  this  study's  domain  is  very 
close  to  this  storm  convergence  area  off  of  the  coast  of  Labrador,  and 
may  cause  the  sharper  than  normal  curvature  of  the  paths  of  the  storms, 
as  suggested  by  the  variability  pattern. 

Figure  5b  is  the  spatial  phase  plot  for  component  three.  The 
change  in  the  relative  phase  angles  over  the  western  part  of  the  plot 
shows  a  progression  in  the  fluctuations  from  north  to  south. 

Meanwhile,  the  variability  area  centered  over  the  coast  of  Labrador  has 
phase  angle  shifts  feeding  in  from  the  Great  Lakes  area  and  from  off 
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the  Atlantic  Coast  around  Newfoundland,  meeting  near  60°N,  60°W,  in 
phase,  and  then  moving  off  to  the  north. 

Over  the  areas  of  low  variability  there  is  a  general  west  to  east 
propagation  of  the  fluctuations  (similar  to  component  one).  These 
fluctuations  are  in  phase  with  the  two  main  areas  of  variability.  This 
is  understandable  since  all  of  these  spatial  patterns  are  connected  to 
the  general  circulation  pattern,  while  there  are  areas  that  exhibit 
greater  spatial  variability  in  pressure  than  others. 

2)  temporal  properties 

Figure  5c  is  the  temporal  amplitude  plot  for  component  three. 

Table  2  shows  that  component  three  was  strongest  in  1976  and  the 
weakest  in  1974.  The  year  of  1976  is  also  the  most  active  year  for 
component  two.  The  fact  that  components  two  and  three  are  most  active 
in  the  same  year  is  not  surprising  since  they  both  seem  to  represent 
polar-outbreak  events,  just  in  different  spatial  areas.  The  Rocky 
Mountains  seem  to  be  the  dividing  line  between  the  two  areas. 

Figure  5d  is  the  temporal  phase  plot  for  component  three.  It 
shows  that  component  three  is  composed  of  a  combination  of  fairly 
rapid,  linear  change  of  phase  events,  and  slower,  sometimes  nonlinear 
change  of  phase  events.  This  fits  in  well  with  the  idea  of  two 
different  types  of  events  making  up  this  component--the  faster  moving 
cyclones  over  the  northeastern  portion  of  the  domain  and  slower  moving 
anticyclones  over  the  western  portion  of  the  domain. 


Figure  5d. 
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3)  comments 

This  component  proves  to  be  the  hardest  to  interpret.  It  appears 
to  be  composed  of  two  spatially  distinct  events  that  cover  a  time 
period  of  6.4  days  to  9.0  days.  The  western  variability  area  is  weaker 
than  the  eastern  one,  but  seems  to  be  more  easily  interpreted  than  the 
eastern  one. 

The  coherence  check  (table  1)  shows  that  there  is  a  significant 
amount  of  coherence  between  component  three  and  the  other  three 
components  being  studied.  We  can  detect  the  presence  of  component  one 
in  component  three  if  we  concentrate  on  the  low  variability  areas  in 
the  spatial  phase  function  of  component  three,  i.e.,  general  west  to 
east  shift.  The  north  to  south  movement  to  the  west  of  the  Rockies  has 
some  movement  over  the  Rockies  as  does  component  two,  so  it's  not  a 
clear  separation  of  the  two  areas. 

e.  component  number  four 

This  component  accounts  for  approximately  7.3  percent  of  the 
variance  in  the  band-passed  data  set.  The  periods  of  the  most 
important  wave  structures  contributing  to  the  spectrum  of  this 
component  were  estimated  to  be  6.9  and  9.0  days. 

1)  spatial  properties 

Figure  6a  is  the  spatial  amplitude  plot  of  component  four.  It  has 
two  main  areas  of  variability,  the  first  being  centered  over  the 
eastern  portion  of  Hudson  Bay  and  the  northern  portion  of  Quebec.  This 
area's  variability  is  almost  equal  to  the  variability  of  the  other  area 
present  and  separated  from  it  by  an  area  of  low  variability.  This 
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second  area  of  maximum  variability  has  a  double  center:  one  located 
near  40°N,  70°W,  and  the  other  near  50°N,  40°W.  These  two  centers  are 
part  of  an  elongated  area  whose  axis  extends  from  Arkansas  through 
North  Carolina  and  Virginia,  and  then  out  along  the  Atlantic  Coast 
through  the  two  centers  mentioned  above. 

This  second  area  off  the  Atlantic  Coast  has  a  close  resemblance  to 
the  "bomb"  deepening  area  as  reported  by  Roebber  (1984)  for  the  years 
1976-1982.  (A  "bomb"  is  a  cyclone  that  deepens  at  least  one  mb  per 
hour  for  24  hours.)  The  works  of  Colucci  (1976)  and  Sanders  and  Gyakum 
(1980)  also  reveal  a  good  match  of  the  offshore  part  of  component  four 
to  the  area  of  maximum  deepening  of  the  intense  extratropical  surface 
cyclones  known  as  "bombs."  Results  of  Zishka  and  Smith  (1980)  show 
this  area  to  be  a  preferred  area  for  "non-bomb"  cyclones  also. 

The  area  of  variability  over  Hudson  Bay  appears  to  be  related  to 
anticyclone  activity.  The  "bombs"  usually  develop  after  the  excursion 
of  cold  Arctic  air  out  of  Canada  down  the  eastern  seaboard.  This 
connection  will  be  discussed. 

Figure  6b  is  the  spatial  phase  plot  of  component  four.  The  phase 
angles  increase  to  the  east  and  to  the  southwest  from  60°N,  102°W.  The 
portion  that  moves  east  corresponds  to  the  area  of  high  variability 
over  Hudson  Bay.  This  area  is  between  120  to  180  degrees  out  of  phase 
from  the  area  off  the  Atlantic  Coast.  These  two  coastal  centers  of 
variability  seem  to  build  in  place  and  die  out  without  much  information 
transferred  between  them. 


2)  temporal  properties 

Figure  6c  is  the  temporal  amplitude  plot  of  component  four.  From 
this  plot  and  table  2  it  can  be  seen  that  component  four  was  most 
active  in  1974  and  the  least  active  in  1981.  In  1981  there  was  a 
trough  at  700  mb  over  the  east  coast  of  North  America  which  forced 
storms  to  track  further  east  than  normal. 

Figure  6d  is  the  temporal  phase  plot  for  component  four.  The  plot 
shows  that  component  four  consists  of  a  large  number  of  slow  phase 
changing  events  and  a  smaller  number  of  rapid  phase  changing  events. 

3)  comments 

It  appears  this  component  picks  up  the  maximum  deepening  events 
(“bombs")  off  the  Atlantic  Coast  and  the  associated  high  pressure 
center  of  the  Arctic  air  which  moves  down  the  eastern  seaboard.  The 
question  to  be  asked  is  why  isn't  the  spatial  phase  plot  showing 
movement  of  the  anticyclone  down  the  east  coast  of  the  US  and  the 
movement  of  the  "bombs"  up  the  coast,  instead  of  what  appears  to  be  in 
place  development? 

The  anticyclone  area  is  cut  off  by  a  low  variability  area  which 
resembles  the  high  variability  area  extension  from  component  two  in 
this  area.  The  shape  of  this  northern  area  in  component  four  therefore 
appears  to  be  artificially  induced  by  the  CPCA  technique. 

The  coherence  check  shows  (table  1)  there  is  coherence  between 
component  one  and  four.  So  component  one  may  have  already  accounted 
for  the  overall  "average"  variability  and  movement  of  cyclones  off  the 
east  coast  of  the  US,  while  component  four  captures  the  most  intense 
deepening  and  movement  stages  of  these  cyclones.  As  noted  by  Sanders 


6.  CONCLUSIONS 


A  data  set  was  constructed  from  eight  years  of  wintertime 
observations  (JAN-MAR)  of  sea  level  pressure  (SLP),  and  then  filtered 
to  eliminate  information  outside  the  2  to  10  day  time  period.  This 
filter  also  augmented  the  original  data  set  with  its  Hilbert  transform. 
The  complex  principal  component  analysis  (CPCA)  technique  was  then 
applied  to  this  augmented  data  set.  This  analysis  technique  is 
designed  to  locate  regions  of  spatially  large  covariability  and  to 
detect  propagating  features  in  the  data  field. 

The  preliminary  interpretation  of  the  first  four  components  of  the 
CPCA  (accounting  for  approximately  63.4  percent  of  the  variance  in  the 
original  data)  are  as  follows: 

1)  Component  one  seems  to  represent  the  variability  pattern  that 
is  produced  by  the  movement  of  anticyclones/cyclones  along  the  mean 
mid-tropospheric  flow.  As  such,  it  seems  to  have  some  of  the  other 
three  components  mixed  in,  and  they  are  all  correlated  as  shown  by  the 
coherence  check  (table  1). 

2)  Component  two  appears  to  show  two  areas  of  variability.  One  is 
along  the  United  States-Canadian  border  and  appears  to  be  caused  by 
anticyclones/cyclones.  The  other  area,  just  to  the  east  of  the 
Rockies,  appears  to  be  linked  to  the  cold  air  outbreak  anticyclones  in 
a  north-south  direction. 

3)  Component  three  seems  to  show  two  distinct  areas  of 
variability.  One  is  possibly  linked  to  the  cold  air  outbreaks  over  and 
to  the  west  of  the  Rockies.  The  other  area  centered  over  the  coast  of 
Labrador  is  harder  to  interpret,  but  appears  to  be  linked  to  cyclone 
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activity  from  the  Great  Lakes  and  off  the  coast  of  the  northeastern 
United  States  and  eastern  Canada. 

4)  Component  four  appears  to  show  the  outbreak  of  Arctic  air  which 
precedes  the  possible  formation  of  a  rapidly  intensifying  cyclone 
called  a  "bomb"  off  the  Atlantic  Coast  of  North  America.  This 
component  also  appears  to  show  the  centers  of  preferred  "bomb" 
formation  areas. 

Limitations  of  the  analysis  technique  appear  to  have  influenced 
the  results  obtained.  While  there  are  rotation  procedures  available  to 
overcome  these  limitations,  they  appear  to  have  their  own  limitations 
and  may  require  the  production  of  a  large  number  of  alternate 
solutions.  More  work  must  be  done  in  this  area  to  eliminate  the 
arbitrary  phase  angle  problem. 

There  is  an  obvious  need  for  further  work  in  the  application  of 
the  CPCA  technique  to  a  SLR  data  set,  and  the  subsequent  interpretation 
of  the  CPCA  results  as  representing  synoptic  scale  disturbances.  A  few 
ideas  for  future  research  would  be  to  1)  expand  the  spatial  area 
covered  to  look  at  the  whole  Northern  Hemisphere;  2)  apply  the 
technique  to  a  larger  number  of  years  of  data  and  partition  the  data 
into  sets  by  season,  and  see  if  the  preferred  areas  shift  with  the 
seasonal  mean  mid-tropospheric,  flow;  3)  apply  the  technique  to  a  set  of 
data  consisting  of  sea  level,  850,  700,  500,  and  300  mb  pressure, 
temperature,  and  height,  and  see  how  the  analysis  of  each  level  and 
variable  fits,  or  doesn't  fit,  the  general  models  of  the  structure  of 
cyclones  and  anticyclones;  and  (4)  reconstruct  an  estimate  of  the 
original  data  set  from  the  "significant"  components  found  by  the  CPCA 
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technique,  and  see  if  the  estimated  data  field  displays  characteristic 
pressure  patterns  and  movement  expected  of  synoptic  scale  disturbances. 

As  is  evident  from  this  study,  the  interpretation  of  CPCA  results 
is  difficult  at  times,  even  with  supporting  information.  Most 
important,  always,  is  that  the  results  make  sense  physically  and  can  be 
found  in  the  original  data  set. 
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ABSTRACT 


ESTIS,  FRANK  L.  Complex  Principal  Component  Analysis  of  Sea  Level 
Pressure  over  Eastern  North  America  and  the  Western  Atlantic  Ocean:  A 
Search  for  Preferred  Disturbance  Activity  Areas.  (Under  the  direction 


of  JERRY  M.  DAVIS.) 

The  purpose  of  this  study  was  to  identify  preferred  areas  of 
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analysis  technique  called  complex  principal  component  analysis.  The 
analysis  technique  was  applied  to  a  sea  level  pressure  data  set 
consisting  of  twice  daily  observations  (00  and  12  GMT)  during  the 
months  of  January,  February,  and  March,  for  the  years  1973-1976  and 
1979-1982.  The  months  of  January  through  March  were  chosen  to 
encompass  the  time  frame  of  the  Genesis  of  Atlantic  Lows  Experiment 
1986.  This  experiment  was  designed  to  study  the  evolution  of  cyclones 
along  the  Atlantic  Coast  of  the  United  States,  because  these  cyclones 
can  become  quite  intense  during  the  winter  months  and  pose  a  threat  to 
ocean  shipping  activities  and  coastal  regions  of  the  United  States  and 
Canada.  The  geographical  area  studied  extended  approximately  from 
30°-12C°W  and  from  25°-60°N,  on  a  polar  stereographic  map  projection. 

The  complex  principal  component  analysis  technique  allows  the 
detection  of  propagating  features  in  the  pressure  data  set.  The 
original  data  set  is  augmented  with  its  own  Hilbert  transform,  which  is 
nearly  equivalent  in  the  time  domain  to  the  quadrature  spectrum  in  the 
frequency  domain.  This  augmented  data  set  thus  contains  information  on 
the  relative  phase  relationships  among  the  propagating  features  of  the 


data  set. 


The  eigenvalues,  complex  eigenvectors  (components),  and  complex 
component  scores  are  extracted  from  the  augmented  data  matrix  and  used 
to  compute  four  measures  that  define  possible  propagating  features  of 
the  data  set.  When  the  component's  values  are  plotted,  at  their 
respective  grid  points,  the  patterns  can  be  seen  that  account  for  the 
largest  amount  of  variance  (first  component),  the  second  largest  amount 
of  variance  (second  component),  and  so  on,  in  the  sea  level  pressure 
data. 

The  variance  in  the  pressure  field  is  felt  to  be  a  direct 
reflection  of  disturbance  activity.  Furthermore,  when  the  data  are 
filtered  to  a  "synoptic  scale"  frequency  window,  the  variance  in  the 
pressure  field  corresponds  closely  with  cyclone  and  anticyclone  storm 
tracks  over  a  time  period  of  2-10  days.  This  spatial  and  temporal 
correspondence  between  synoptic  scale  disturbances  and  variance  of 
pressure  allows  the  use  of  the  complex  principal  component  analysis 
technique  on  the  filtered  sea  level  pressure  data  set  to  objectively 
find  preferred  areas  of  disturbance  activity. 

The  results  of  the  analysis  produced  four  components  that  could  be 
given  a  preliminary  interpretation  from  a  physical  sense.  The 
components  (patterns)  resembled  areas  of  activity  as  discovered  in 
previous  studies  using  other  analysis  techniques.  Component  one 
resembled  the  variability  pattern  produced  by  the  movement  of 
anticyclones/cyclones  guided  by  the  mean  mid-tropospheric  flow. 
Component  two  represented  the  cold  air  outbreak  anticyclones  to  the 
east  of  the  Rockies,  and  anticyclones/cyclones  which  move  along  the 
United  States-Canadian  border.  Component  three  is  linked  to  the  cold 


air  outbreak  anticyclones  to  the  west  of  the  Rockies,  and  to  cyclone 
activity  over  the  Great  Lakes  and  off  the  coast  of  the  northeastern 
United  States  and  eastern  Canada.  Component  four  represents  the 
formation  of  the  "bomb"  cyclone  off  the  Atlantic  Coast  of  North 
Carolina,  and  the  preceding  Arctic  air  outbreak  anticyclone.  The 


